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Abstract 

Often considered in numerical simulations related to the control of quantum systems, 
the so-called monotonic schemes have not been so far much studied from the functional 
analysis point of view. Yet, these procedures provide an efficient constructive method 
for solving a certain class of optimal control problems. This paper aims both at ex- 
tending the results already available about these algorithms in the finite dimensional 
case (i.e., the time-discretized case) and at completing those of the continuous case. 
This paper starts with some results about the regularity of a functional related to a 
wide class of model in quantum chemistry. Those enable us to extend an inequality due 
to Lojasiewicz to the infinite dimensional case. Finally, some inequalities proving the 
Cauchy character of the monotonic sequence are obtained, followed by an estimation of 
the rate of convergence. 
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1 Introduction 

Following the increasing interest of the chemists community for optimal control of quantum 
systems [21, 27] and the successful laboratory demonstration of control over molecular phe- 
nomena (see, e.g., [1, 7, 32] and more recently [12, 31]), some mathematical studies of the 
models involved in this topic have been carried out, see e.g. [6, 16]. In this way, it has 
been proved in recent papers [4, 9] that a wide class of optimization problems considered by 
chemists are well posed. Yet, these proofs are not constructive and consequently do not give 
rise to concrete numerical methods to approximate their solutions. 

On the other hand, at numerical simulation level [8, 22], various kind of procedures 
exist and show a good efficiency. Among them, the so-called monotonic algorithms have 
demonstrated their efficiency on several problems. In a recent paper, a study of the time- 
discretized algorithms [24] have been presented and first functional analysis results have been 
obtained about the continuous case [14, 25]. 

The aim of this paper is to complete these works by providing general proofs of conver- 
gence of the optimizing sequences. Consequently, we obtain a constructive method, indepen- 
dent of time or space discretization to compute critical points (and sometimes extrema, see 
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Remark 3) of the cost functional under consideration. 

Let us briefly present the monotonic schemes in the simple case of ordinary differen- 
tial equations (ODE). Let A, B,C be three square matrices in .M n (R), C being symmetric 
positive, a > and T > 0. Consider the optimal control problem corresponding to the 
maximization of the functional J defined by: 



T 

2/ 



J(v) = y(T) ■ Cy(T) - a v\t)dt, 



f 

— a I v 
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where "•" denotes the usual scalar product of W\ Here, the state y : [0, T) — > R™ and the 
control v : [0, T] -» R are linked by the ODE: 

y'(t) = (A + v(t)B)y(t), Vte(0,T) 
2/(0) = 2/o 

the initial condition yo being fixed. 

Given two controls v and u and the corresponding states y and y, we first note that: 
J(v)-J(v) = (y(T)-y(T))-C(y(T)-y(T))+2(y(T)-y(T))-Cy(T) 

-a (v(t) -v(i))(v(i) +v(t))dt. 
Jo 

We then introduce an auxiliary function z : [0, T] — > R™ associated to y and d by 

*'(t) = -(A*+t;(t)S*)a;(t), 
= Cy(T) 

where A* and B* are the transposed matrices of A and B. 

Focusing on the second term of the right hand side of this equation, we get: 



(y(T) - y(T)) ■ Cy(T) = f (v(t) - v(t))By(t) ■ z(t)dt. 

Jo 

Thus, we finally obtain: 

J(?0-J(t;)=( 2 7(T)-y(T)).C(y(T)-y(T))+a^ T (^)-^W) ' «(*)"«(*) - < t )) dL 

A simple way to guarantee that v gives a better cost functional value than v, is to impose 
that: 

(v(t) - v(t)) (^By(t) ■ z(t) - v(t) - v(t)j > 0. (f ) 

Following this approach, the sequence (v k )keN defined iteratively by the implicit equation 
y k+i _ ±By k+1 (t)-z k (t), where y k+1 andz fc correspond to v k+1 andw fc respectively, optimizes 
J monotonically since 

J(v k+1 ) - J(v k ) = (y k+1 (T) - y k (T)) ■ C (y k+1 (T) - y k (T)) + a f (v k+1 (t) - v k (t)) 2 dt > 0. 

Jo 

In this article, we prove the convergence of generalizations of this algorithm towards a critical 
point of J in the case of the Schrodinger partial differential equation: 

idtip(x, t)-[H - n{x)e{t)]^{x, t) = 0. 

This equation governs the evolution of a quantum system, described by its wave function 
■0, that interacts with a laser pulse of amplitude e, the control variable. The factor fj, is 
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the dipole moment operator of the system. In what follows, H = —A + V where A is the 
Laplacian operator and V = V{x) the electrostatic potential in which the system evolves. 
We refer to [22] for more details about models involved in quantum control. 

The paper is organized as follows: we start in Section 2 with some necessary results about 
the linear and nonlinear Schrodinger equations involved in the problem we are considering. 
We then present the optimization problem in Section 3, and claim some regularity results 
about the corresponding cost functional in Section 4. We introduce in Section 5 an important 
tool for proving the convergence of the sequence, namely the Lojasiewicz inequality and some 
of its generalizations. The definition of the monotonically optimizing sequence is given in 
Section 6 where some useful properties are also claimed. The convergence of the sequences 
is proved in Section 7 and a first result about their rate of convergence follows in the last 
section. 

Throughout this paper, T is a positive real number representing the time of control of a 
physico-chemical process. We denote by L 2 and L°° the spaces L 2 (R 3 ,C) and L°°(R 3 ,C), 
W p -°°(R 3 , M) with pe[l, +oo) by W p '°°, the Sobolev space H 2 {R 3 , C) by H 2 and U>(0, T; X), 
with p e [l,+oo) denotes the usual Lebesgue space taking its values in a Banach space X. 
We also use the notation ( . | . | . ) and ( . , . ) defined by: 

(f\A\g)= [ J(x)Ag(x)dx, (f,g)= [ J(x)g(x)dx, 

where / and g are in L 2 and A is an operator on L 2 . To simplify our notation, the space 
variable x will often be omitted. Finally, for h e L p (0,T; X), p e]l,oo], we recall that 
II^IIlp(o.t ; x) = \\t !— » II M*) II -HI lp (o t) ■ Finally? we denote by Im(z) and Re(z) the imaginary 
and the real part of a complex number z. 



2 Preliminary existence results 

The sequences we study in this paper are defined through iterative resolutions of Schrodinger 
equations. Before introducing the relevant framework of our study, we present here some 
necessary preliminary existence and regularity results concerning these equations. The first 
one will correspond later to the initialization step in the definition of the sequences. This 
lemma is a corollary of a general result on time dependent hamiltonians (see [23], p285, 
Theorem X.70) but for the sake of clarity, we present here an approach using other techniques 
also useful in the proof of the next lemmas. 

Lemma 1 Let fi and V belong to W 2 -°° and let H = -A + V. If e G L 2 {0,T) and ip € H 2 , 
the equation 

f id t iP(x, t) - [H(x) - n(x)e(t)]iP(x, t) = . . 

\ ^(x,0)=Mx) 1 ' 

has a unique solution ip e L°°(0, T; H 2 ) n W llO °(0, T; L 2 ). Moreover: 

Vte[0,T], W{t)\\ L * = \\ML>- (3) 

Proof: One can also read a similar proof in [3] but we give here some details. It is well 
known (sec [II] for instance) that for any T > and u a € H 2 , the Schrodinger equation 

f id t u(x, t) + Au(x, t) = 0, iel, te[0,T] 
\ u(x, 0) = uo(x), x e R 

has a unique solution u(t) = S(t)u such that u G C([0,T];H 2 ) n C 1 ([0, T]; L 2 ), where 
(S(t))t£R denotes the free Schrodinger semi-group e ltA . Moreover, for all t € [0,T] we have 
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\\u(t)\\ H 2 = \\S(t)u Q \\ H 2 = ||uo||jj2. (4) 

Let A > be a given positive number which will be fixed hereafter and denote by Y the 
space C([0, T]; H 2 ) endowed with the norm \\i/j\\y = sup tg [ T i e _At ||'i/'(i)ll-H' 2 - The solution of 
equation (2) is obtained equivalently as a solution to the integral equation 



= S(t)i)Q+i( S{t- s)W(s)^(s)ds 
Jo 



where W(x, t) = —V(x) + fi(x)e(t) for all t e [0, T], x G R 3 . We are going to show that this 
equation has a unique solution in Y, by proving that operator $ defined by 

$(1>)(t) = S(t)ifo + i f S{t- s)W{s)ip{s)ds 
Jo 

has a unique fixed point in a closed ball Br = {tp £ Y ; \\iP\\y < R} for suitable R. 

If V € B R , then ||^(s)||ff= < e As ||?M|y < Re Xs and since e £ 2 (0, T; VK 2 <°°), we can set 
(0 > such that || VF||l2( ,t;W 2 .°°) < P- Using estimate (4) and Cauchy-Schwarz inequality we 
obtain 

i 

t / rt 

2\s 



\Mm)\\ H i < \\Mh* + y o \\w(s)Tp(s)\\ H2 ds < \\ik\\H*+pR{j ^ s ds 

R 

It follows that if R > is large enough so that H^oIIh 2 < ^ and if we choose A > 2p 2 , then 

i_ 

||*W0lk< sup e- xt \\M\H*+pR{ fo 2X ^dX <| + < fl. 

This means that $ maps B/j into itself. Then, for tjji , ip 2 € £?_r it is clear that 

/e 2At - f 

||(*(V>i)-^2))WII^ < y o ||W( a )(^i-^)(«)||ffid*<p||Vi-^l|y 

and since A has been appropriately chosen, this proves that <I> is a strict contraction from Br 
into itself as 

- < P llV-i - iMv sup f 1 " 6 / < -^=||Vi - H\y < h\1>i - My 

te[o,T] V 2A / V2A 2 

and therefore $ has a unique fixed point, yielding the solution of equation (2) in L°°(0, T; H 2 ). 
One can notice that uniqueness is not only true in Br but also easily proved using the norm 

in L°°(0,T; L 2 ). Moreover, calculating Ira J (2).ip(x) dx, one can prove the conservation of 
the L 2 -norm (3) and finally, using equation (2), it is easy to obtain that tp £ W 1,oo (0, T; L 2 ). 



We will also have recourse to a similar lemma, dealing with equation (2) with a non zero 
source term. 

Lemma 2 Let H, /j,, e be as above and ip £ L°°(0, T; H 2 ). Given e' € L 2 (0, T), the equation: 

/ id t tp'(x, t) - [H{x) - ii{x)e{t)]i)'{x, t) = -^{x)e' {t)^{t, x) 

1 V'(x,0)=0 K ' 
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has a unique solution ip' e L°°(0, T; H 2 ) n W 1 ' 00 ^, T; L 2 ). Moreover i/ie following estimate 
holds: 

\\^'\\l=°(0,T;L 2 ) < 2||^||lHI £ 1U 1 (0,T)||V'I|l~(UT;L 2 )- (6) 

Proof: The key point to prove the existence of a solution for (5) is to underline the fact 
that the source term f(x, t) = /j,(x)e'(t)ip(t, x) of this linear Schrodinger equation belongs to 
L 2 (0,T;H 2 ). It is then very classical to get from the Lemma 1 the existence and uniqueness 
of a solution V' to equation (5) in L°°(0, T; H 2 )r\W 1,o °(0, T;L 2 ). Consider now <p E C([0,T\) 
defined on [0,T] by ip(t) = \\4>' (t)\\ 2 L2 . We have: 

j t <p(t) = 2Re (yj'(t), g ~^ £( *V (*) - = '^'(t) Im (V' • (7) 

Moreover, there exists to such that: <p(to) = sup 4g [ r i { || ip'(t)\\ 2 L2 } • We suppose e' ^ and 
ip 7^ 0, so that <p(to) 7^ by uniqueness of the solution of (5). Since ip'(x, 0) = for all x £ K 3 , 

/•to 

the integration of (7) between and to yields ip(to) = / — 2e'(t)Im{tp'(t)\fi\ip(t))dt , then: 

Jo 

<p(to) = \mto)\\h < ml-wwwv [ T 2\e'(t)\\\m\\Lidt. 

Jo 

Since ||^'(*)||l= < ||V>'(to)||z,2 for all t e [0,T], we obtain 

||^(t)IU» < ||V>'(*o)IU* <2||m||l-||VIUoo( 0;W [ T \e'(t)\dt, 

Jo 

what ends the proof of estimate (6). 

■ 

Finally, we claim a last result that will be useful to tackle the problems related to a nonlinear 
Schrodinger equation we encounter in this study. Actually the nonlinearity we consider here 
is the one that appears naturally in the adjoint system from a quadratic cost functional (as 
J is in (14)), even when the state equation is linear. 

Lemma 3 Let H, fj,, e and ipo defined as above. Given \ € L°°(0, T; H 2 ), the nonlinear 
Schrodinger equation: 

f idrt(x,t)-[H(x)-v(x)e(t)+Irn( X {t)\ f im)) f i(x)}iP(x,t)=0 ( 
\ #c,0)=Vo(z) {0> 

has a unique solution ip G L°°(0, T; H 2 ) n W llOO (0, T; L 2 ). 

Proof: - First Step - 

Let u and \ € H 2 , we denote the nonlinear term by F(u) = Im{x\lA u )l JLU an( i we can P r ove 
that one has the following estimates: 3C = C{\, fi) > such that 

Vu,veL 2 , \\F(u)-F(v)\\ L 2<C(\\u\\ L 2 + \\v\\ L 2)\\u-v\\ L 2 (9) 
Vu,veH 2 , \\F(u) - F(v)\\ H 2 < C(\\u\\ L 2 + \\v\\ H 2)\\u-v\\ H 2 (10) 

\\F( U )\\ H 2 < C\\U\\ L 2\\ U \\ H 2 (11) 

Indeed 

\\F(u) - F(v)\\ L 2 < \\Im(x\n\u)nu- Im(x\fJ-\v)nv\\ L2 

< \\Im(x\v\u)n(u- v)\\ L2 +\\Im{x\fi\(u-v))fj,v\\ L2 

< IHIi,°°llxlU 2 (IMU 2 + \\v\\ L2 )\\u- v\\ L 2 
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which proves (9). Now, we have to establish (10) and (11). First of all we have 

\\F(u) - F(v)\\ 2 H2 = \\F{u) F{v)\\\ 2 + \\AF(u) AF(v)\\ 2 L2 . 
The first term of the right hand side is conveniently bounded in (9). Moreover 

||AF(«) - AF(v)\\ L 2 < \\Im{ X \»\u - v)A(»(u - v))\\ L2 + \\Im( x \^\v)A^v)\\ L2 

< ||m||^2,oo||xIU 2 (II u IIl2 + \\v\\ H 2)\\u- v\\ H 2 

< C(\\u\\ L 2 + \\v\\ H 2)\\u- V\\ H 2. 

Then, F is locally lipschitz in H 2 . Therefore, taking v = 0, we also get (11). 
- Second Step - 

The proof of a local-in-time result is based again on a fixed point theorem. We begin by 
fixing an arbitrary time T > and considering r e]0,T]. We also consider the functional 



£ : V — > U( . , 0)^o -i f'u(., 8)F^(s)) ds, 
Jo 



where {U(t, s), s,t <E [0,T]} is the propagator associated with the operator H — [is and 
induced by Lemma 1 (such that U(t, s) E £(H 2 ) - for details, see [5]), and the set 

B = {v e L°°(0,t; H 2 ), ||V|Uoo (0 , r; ^) < 2M||V>o||W- 
where M satisfies Vw e H 2 , \\U(t, s)v\\ H 2 < M\\v\\ H 2. 

If t > is small enough, the functional £ maps B into itself and is a strict contraction in 
the Banach space L°°(0, r; H 2 ). Indeed, on the one hand, from estimate (11), if ip € B, we 
have for all t G [0, r]: 



mm)\\ 



H 2 



< 



-i f U(t,s)F(^(s))ds 
Jo 



H 2 



< M\\Mh 2 +tM\\F WWlooq^h*) 

< M\\^\\ H 2 +rCM||V'|U°o(0,r;I, 2 )||V'IU°°(0,r;H 2 ) 

< M||Vo||jf2 +A T CM 3 \\i; Q \\ 2 H2 . 

Then, if we choose r such that 4tCM 2 \\iPo\\ H 2 < 1 we obtain ||^(V')l|L° o (o.T;_ff 2 ) < 2M||^o||j? 2 
and belongs to B. On the other hand, if ipi and ^2 € B, then for all t in [0, r] we have, 



\mi)(t)-a^)(t)\\ 



H 2 



H 2 



U(t,s) (F(M*)) ~ F(M*))) ds 

< CM (||'0l||i° o (O.-r; J L 2 ) + II^2||l=°(0.t;_H" 2 )) / || V>1 («) ~ ^2 («) \\ H 2 ds 

Jo 

< 4tCM 2 \\iI>o\\ H 2 - 1p2\\L™(0,T;H 2 ) 

with 4tCA^ 2 1 1 -0o 1 1 jf 2 < I- Therefore, from a usual fixed point theorem, we can deduce 
existence and uniqueness in the set B, then in L°°(0, r; H 2 ), for r > small enough, of the 
solution of equation 



-/ 

Jo 



U(t,s)F{ip(s))ds 



(12) 



which is in fact equivalent to equation (8). Moreover, using (8), it is easy to prove that d t ip 
belongs to i°°(0,r;i 2 ). 



G 



The last point is then to prove the uniqueness of the solution u of (8) in the space 
L oo (0,T;i7 2 )nVK 1 ' oo (0,T;L 2 ). Let Va and ^2 be two solutions of (8) and w equal to ipi - V>2- 
Then w(0) = and 



id t w - [H(x) - fj,(x)e(t)]w = F{tp 2 ) - 



(13) 



Calculating Im 



(13).w(x) dx and using ( !)) we obtain j^dlHIi 2 ) — ^ll^lli 2 an< ^ unique- 



ness follows by Gronwall lemma. Hence the proof of uniqueness, existence and regularity of 
the solution of equation (8) in R x \0,t] for any time r < — — r „.. , , — . 

- Third Step - 

Now, the goal is to obtain an a priori estimate of the solution in W 1,oo (0, T; L 2 )nL°°(0, T; H 2 ) 
for any arbitrary time T, in order to prove that the local solution we obtained previously 
exists globally because we have a uniform bound on the norm H^WII h 2 + ||^tV'(*)IU 2 - 

Actually, since equation (8) is equivalent to the integral equation (12) and since it is easy 
to prove the conservation of the L 2 -norm of the solution, we have, 



\M)\\h' 



< 



< 



\\U(t,0)i> \\ H2 + [ t U(t,s)F(,p(s))ds 
Jo 

M\\Mm+MC f U{s)\\ L ,U{s)\\ 
Jo 

t (l + J*\ms)\\ H * ds 



H 2 

ds 



< Co 



where Cq.t > is a generic constant depending on the time T, on /i, \ and on H^oll^ 2 - We 
finally obtain from Gronwall lemma and from equation (8), that ||V>(*)l|i/ 2 + ll ( 9t' ! / ; (0IU 2 ^ Co,t 
for all t £ [0, T]. Hence the proof of Lemma 3. 



3 Optimization problem 

Let us now present the optimization problem we are dealing with in this paper. Let O be a 
positive symmetric bounded operator on H 2 and a and T two positive real numbers. Given 
ipo £ H 2 , we consider the cost functional J defined on L 2 (0,T) by: 

J(e) = (V(T)|0|V>CT)) - " [ T e 2 (t)dt, (14) 

Jo 

where tp is the solution of (2). In all the sequel we suppose that HV'oIIl 2 = 1. The existence 
of a minimizer for similar cost functionals (with the opposite sign) has been obtained in [3] , 
[4] and [9] and follows from the construction of a minimizing sequence and a compactness 
lemma (Aubin's lemma). Here, the point is to maximize the functional J and as usual, at 
the maximum of J, the Euler-Lagrange critical point equations are satisfied. A standard way 
to write these equations is to use a Lagrange multiplier i) usually called adjoint state. 
The following critical point equations are thus obtained, for i£K 3 and t £ (0, T): 

id t tp(x, t) - [H(x) - fj,(x)e(i)]i>(x, t) = 0, , . 

( id t x{x, t) - [H(x) - fi(x)e(t)]x(x, t) = 0, . . 

\ X (x,T) = Oi,(x,T), 
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ae(t) + Im(4>(t)\v\x(t)) =0. 

The existence of \ € £°°(0, T; iJ 2 ) results from an adaptation of Lemma 1, as for ip{T) e H 2 
since equation (15) is actually equation (2). In what follows, we also consider the linearized 
equation of (16): 

id tX '{x, t) - [H(x) - n(x)e(t)] X '(x, t) = -/i(x)e'(t) X (t, x) 
X '(x,T) = 0^'(T), {U) 

where e' € L 2 (0,T) and ip' is the solution of (5), corresponding to the solution ip of (15). 
The existence of x' € L°°(0,T; H 2 ) follows from Lemma 2. The analysis done in the proof 
of estimate (6) gives in this case: 

||x'(i)|| L2 < 2||/i|| i oo|| e '|| i i ( o, r) ||x||L»o(o,r;i») + \\X'(T)\\ L 2. (18) 

Since X (T) = Oip(T), x'(T) = Oip'(T) and from (6) and the conservation of the L 2 -norm, 
we obtain 

||X'IU~(0,T;L») < 2||/i||£oo||e'|| i i (0ir) ||x|U<»(0,r;£>) +2||O||*|| / u||l~||£'|Ui ( 0,T)||V'IU-(0,T;L2) 

< 4|| M || i =»||O||*||e / || L i (0 ,T), 

(19) 

where ||0||* denotes the operator norm of O on L 2 . 

4 Properties of the functional J 

We begin with some properties about the regularity of the cost functional J. 
4.1 Gradient of J 

We start with some first order properties. As often, the use of the adjoint state x allows us 
to simplify the computation of the derivative of J. This result is the purpose of the next 
lemma. 

Lemma 4 The cost functional J is differentiable on L 2 (0,T) and its gradient can be ex- 
pressed by 

(VJ(e), e') = -2 [ (ae(t) + Im{ X (t)\^(t)))e' (t)dt, (20) 
Jo 

where (•, •) is the usual inner product on L 2 (0, T) and ip and x are the solutions of (15) and 
(16). 

Proof: We only give here a sketch of the proof. The details can be found in reference 
[3] for a slightly different cost functional. The main point is to prove the differentiability of 
the functional <f) ■ £ € L 2 (0,T) 4>{T), where ip is the solution of equation (15). Actually, 
one can prove that the solution ip' of (5) is such that D<j>(e)[e'] = ip'{T). Therefore, since 

J( e ) = (ip(T)\0\ip(T)) -a [ e 2 {t)dt, we have 

(VJ(e).e') - 2 Re(ip'{T)\0\ip{T)) -2a [ e(t)e'(t)dt. 

Jo 

To end the proof of (20), we consider the solution x 01 the adjoint state equation (16) and 
we multiply equation (5) by x (the complex conjugate of x), integrate on M x [0, T] and take 
the imaginary part. We obtain: 



Ira I f {idtip 1 -\H- Lie]iP')x = Im f f tie' 

JO JR JO JR 



ipX- 
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After an integration by parts and since ib'(0) = 0, we get 

r T r f i-T 



Im f f idtx^p' -Im [ tp'(T) i X (T) - Im f f [H - iis\xi>' = Im [ f /iety*. 
jo Jr Jr jo Jr. Jo Jr 

Since \ satisfies equation (f6), we then obtain 

Re{ib'(T)\0\^{T))=Re f ^(T)O^T) = -Im f [ »s'ib X =- f /m<x(t)|/i|^(*)>e '(t)dt 
Jr Jo Jr Jo 



what ends the proof of the lemma. 



In what follows, we denote by VJ(e) the function t — 2(ae(t) + Im(x{t)\fi\ip(t)}) and 
by Cj the set of the critical points of J, i.e., 

Cj = {e £ L 2 (0,T), Vie[0,T], ae(t) + Im{x{t)\n\m) = 0}- (21) 

Note that, thanks to the results of the section 2, we have Cj C L°°(0, T) since for all e £ Cj, 

||e|U°»(0,T) < ^||(xM^)l|l/»(0,T) < C\\h\\ L ^\\x\\l^(0,T:L^H\\l^(0,T:L^- 

Remark 1 : Note also that for a > 6T||/x||| i0 o||0||*, the set Cj is reduced to one point. 
Indeed, suppose that Cj contains two distinct points S\ and £2, we then have, for t £ (0,T); 

a(e 2 (t) - ei(i)) + Im( X2 (t) - XiWI^I^W) + Im{ X i(t)\n\Mt) ~ M*)) = «, 

where ipi, ip2 (resp. Xi> Xi) are the solutions of (15) (resp. (16) ^ corresponding to E\ and £2 
respectively. Using estimates (6) with ip = ipi, ip' = 1P2 — tpi, £ — £2 and e' = £2 — £\ and 
(18) with x = Xi? x' = Xi — Xi> i>' =tp2 — ipi, £ = £2 and s' = s 2 — e\, we obtain 

ot\\e2 - ei\\ Ll{0X ) < 6T||/i||io.||0||.||e2 - eilUno.r), 
which leads to a < 6T||0||*||/u||f /00 , and the result follows. 

In order to prove the compactness of Cj, we introduce an important property of the 
application e(t) ip(x,t), firstly presented in a more general setting by J. M. Ball, J. E. 
Marsden and M. Slemrod in [2]. In our context, this result can be stated as follows. 

Lemma 5 Assume that e £ L 1 (0,T), \i : X — > X is a bounded operator and that H generates 
a C° -semigroup of bounded linear operators on some Banach space X. For x £ M 3 and 
t £ (0, T), we denote by ip(x,t) the solution of 

id t ip — [H — /ue]V> = 0, 
V>(0) = tp Q £X. 

Then, e 1— > ip is a compact mapping in the sense that for any weakly converging sequence 
(en)nGN to £ in -L (0, T), (ip n ) n eN converges strongly to ip in C([0, T];X). 

The precise proof of this result derives directly from [2] (Theorem 3.6, p580), see also [10] 
and [26]. It allows us to obtain the following lemma. 

Lemma 6 For fj, £ W 2 >°° , Cj is compact in L°°(0,T). 
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Proof: Consider a bounded sequence (e™) n eN of Cj. By definition, for all n £ N, 

e n £ L 2 (0,T) and £™(i) = (x"(OI/ i IV' n (0) where ^™ and \ n are the corresponding 

a 

solutions of (15) and (16). It is also possible to extract a weakly convergent sub-sequence 
in L 2 (0,T), still denoted (e n ) n eN- From Lemma 1, one knows that the Hamiltonian H = 
—A + V with V £ W 2 - 00 generates a C°-semigroup of bounded linear operators on the Banach 
space X = H 2 . Therefore, with \x £ W 2 -°° the conditions of Lemma 5 are fulfilled and we 
obtain the strong convergences ip n "^-J 30 ^ an d ^™ n-+_oo ^ c([0,T]; H 2 ). Thus, for 

all t £ (0,T), / iP n (t) ii X n (t)dx / iAf) nx(t)dx. The sequence (e„(i)) neN then 

strongly converges in L°°(0,T) and the result follows. 



4.2 Analyticity of J 

The implicit formulation of the derivative can be iteratively carried on in order to prove the 
analyticity of J. 

Lemma 7 Let ip be the solution of (15) corresponding to e. The functional 

i?:L 2 (0,T) -» L oo (0,T;iJ 2 )n^ 1 ' oo (0,T;i 2 ) 

£ i ► ip , 

is analytic. 



Proof: Let e,e' e £ 2 (0,T) be such that ||£'||z,i(o,T) < 



4||/i|U° 



and the sequence 



(^)<€N € (i°°(0, T; ff 2 )) N defined recursively by ip° = d(e) and for I > 0: 
ifit^z, t) - [fT - fi(x)e(t)]ip e (x, t) = -fiixyit)^- 1 ^, t) 

ip e (x,o) = o. 



(22) 



The existence of ip 1 is consequence of Lemma 2. Thanks to (6) applied with ip — ip 1 1 and 
ip' = ip e , one has for £ > 1 and i e [0, T]: 

||^(*)IU" < 2||mII^I|£'IIl 1( o,t)||^- 1 ||l~(o,T;L 2 ) < 2 e MU\e'\\{ H0 ,Ty (23) 
Given AT > 0, we obtain by summing (22) from I = to N: 

/ n \ / n \ 



AT 



\£=0 



U=0 



(24) 



E^(^0) = ^ o (x). 



On the other hand, one has: 

id t d{e + e') - [if - n{x){e{t) + e'(t))]0(e + e') = 0. 



(25) 



Subtracting (24) and (25) and using estimates (6) with ip — —ip N , ip' — Y^t=o i>^( x > *) — ^(£ + 
e'), £ = e + e' and (23), we get: 



N 



^e + e')(t)-J2^(t) 



1=0 



<2 N MUe'\\%, T) <2- N 



L 2 
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and the functional d reads now: tf(e + £ ') = V in L 2 (0, T). Since e' h-> ip e is ^-linear, the 

e=o 

theorem follows. 

■ 

The next lemma follows immediately from this result. 
Lemma 8 The cost functional J is analytic on L 2 (0,T). 

4.3 About the Hessian operator of J 

Let us now investigate some properties of the second order derivative of J. Though we 
express it as an implicit function of its argument e, some results can be obtained from the 
next lemma. 

Lemma 9 Let ip and \ be the solutions of (15) and (16). The functional 7 : e 1— > Im(x\^\ip) 
is differ entiable on L 2 (0,T) and one has: 

D 7 (e)[e'} = Im{ X '\m + Mx W>> (26) 

where ip' and x' are the solutions of (5) and (17). Moreover, for all e G L 2 (0,T), Dj(e) is 
compact on L 2 (0,T). 

Proof: Let e G L 2 (0,T) and ip and x the corresponding solutions of (15) and (16). As 
in the proof of Lemma 4, the key-point is the differentiability of the functional defined in 
Lemma 7 on L 2 (0,T). Actually, D-d(e)[e'] = ip', where ip' is the solution of (5). The main 
explanations can be read in [4]. Repeating this argument for e n Xi we obtain that 7 is 
diffcrcntiable and we get (26). 

Let us now prove the compactness of this operator. Let (£'™)neN be a bounded sequence 
in L 2 (0,T) and let (tp' n ) ne ^ and (x' n )neN be the corresponding solutions of (5) and (17). 
As ip' n € L°°(0,T;H 2 ) n W^°°(0,T;L 2 ) (see the proof of Lemma 2), we have that ip /n S 
C([0,T];L 2 ) and d t ip' n G L 2 {0,T;L 2 ). By means of the continuity of: 

L 2 (0,T) C([0,T];L 2 ) and L 2 (0,T) L 2 (0,T;L 2 ) 

e' ^ ip' s' ^ d t tp' { > 

there exist ip' 00 such that, up to extraction, ip' n ip' 00 G L 2 (0,T;L 2 ) and 

dt^^dt^ eL 2 (0,T;L 2 ). (28) 

Since ip' n (0) = 0, we have ip' n {t) = [ d t ip' n {s)ds and (28) implies that for all t G [0,T], 

Jo 

(||V' / ™(t)||L2) neN is uniformly bounded. Moreover, for all t,t' G [0,T],i < t' , we have: 

||^"(t') - lT (t)|| ia < jf H^WIU'ds < VF^t||9^'"|| i2 (0,T;^)- 

Combining this with (28), we find that (V>'™)neN is an equicontinuous sequence in C([0, T], L 2 ). 
We conclude by applying Ascoli's theorem to the family {Im(x\n\ip' n ) , n G N} of the space 
C([0, T]). Similar arguments apply for {Im(x' n \f-i\ip) , n G N}, and the results follows. 



Thanks to the previous lemma, J is twice differentiable and its Hessian operator reads: 

Hj{e) : s' ^ -2(ae' + Dj(e)[e']). 

In the sequel, a criterion ensuring that the Hessian operator of J is invertible will be useful. 
The next lemma provides it. 
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Lemma 10 Suppose that: a > 6T||/i|||oo||0||*. Then the operator Hj(e) is invertible on 
/. J ;0. / :.. 

Proof: We keep the notation of Lemma 9. The Cauchy-Schwarz inequality, combined 
with (6) and (19) yields: 

p7(e)[e']IU~(0,T) < |MMUx'IU~(o, W + ||O||*||V'|U~(0,T;^)) 
< 6>/r||A*||ico||O||.||e , || i2 ( ,D. 
Finally thanks to the assumption of the lemma, one has 

sup (U\D 7 (e)[e'}\\ L2(0tT) ) <1, 

which implies that / + —Dj(e) is invertible and the result follows. 

a 



5 Lojasiewicz inequality for the cost functional J 

Several convergence results of dynamical systems have been proved thanks to the Lojasiewicz 
inequality recalled here. In order to tackle the problem of the convergence of the optimizing 
sequence presented in the next section, we have to extend this inequality to the case of a 
compact set in an infinite dimensional space. The basic result considered in this section is 
the following (cf [17, 18]): 

Theorem 1 Let N be an integer and T : R N — > R be an analytic function in a neighborhood 
of a point a G R N . Then there exists o > and 9 G]0, \] such that 

\fx e R N , \\x-a\\ < cr, ||vr(x)|| > |r(x) -r(a)| 1 - £ ', (29) 

where ||.|| is a given norm on R N . 

The real number 9 is a Lojasiewicz exponent of a. Following the work [15] of M. A. Jen- 
doubi (which simplifies the theorem of Lojasiewicz-Simon [28]), the latter theorem can be 
generalized to the case of infinite dimension. 

Lemma 11 Given e G L 2 (0, T), there exists a' > 0, k > and 9' G]0, \\ such that: 

Ve' e L 2 (0,T), ||e' - e|| i2(0 ,r) < a*, || V J{e')\\ L ^,T) > k\J{s') - Jie)] 1 - ' . 

We give the proof of this lemma in the appendix. A more precise result can be obtained 
if the Hessian operator under consideration is invertible at point a (see e.g, [13]). Indeed, 
one can then show that 1/2 is a Lojasiewicz exponent of a. We will use this improvement 
in Section 8 since Lemma 10 provides actually an expected sufficient condition. The next 
lemma is a global version of the previous one. 

Lemma 12 Let Cj be a connected component of Cj in L 2 (0,T). We denote by I the value 
of J(e) for all e e Cj and we set J(e) — I — J{e). There exist a > 0, k > and 9 G]0, |] 
such that: 

VeeL 2 (0,T), d 2 (e,Cj) <n, || VJ(e)|| L 2 (0;T) > k\ J{e)\ l ~^ , (30) 
where di is the distance associated to the L 2 (0,T)-norm. 
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Proof: Lemma 11 ensures that for each point a in Cj there exist three real numbers o~ a , 
9 a and n a such that: 

Ve e R N , ||e - o|| L 2 (0 ,D < *a ||V J(e)|Uw) > n a \ Jie)^ . 

The compactness oiCj, guaranteed by Lemma 6, allows us to extract from |-B(a, a G Cjj 
a finite family A = {B(a,i, ~2~)} ieF , where F is a finite set of indexes, such that Cj C A. 

We then define a, k and 9 e]0, 1/2] as the respective lower bounds of {~^ L } ieF , { K ai} ie F 
and {9 ai } i£F and the result follows. 



6 Optimizing sequence 

We have now gathered all the necessary results to present and analyze the optimizing se- 
quence. 

6.1 Definition of the sequence 

Following the approach sketched in the introduction, Y. Maday and G. Turinici have defined 
an optimizing sequence (e k )ken for the cost functional J as follows [20] : 

Consider (5,7]) e]0, 2[x]0, 2[, e° G L°°(0,T), e° € L°°(0,T), V>° and x° the corresponding 
solutions of (15) and (16) according to Lemma 1. The functions e k and e k are computed by 
solving iteratively: 



id t ^ k {x,t) = (H(x)-e k (t)»(x))iP k (x,t) 
i> k (x,0)=Mx) 



e k (t) = (1 - 5)e k -\t) -Im( X k -\t)\n\i, k {t)) (32) 

a 



id tX k (x,t) = (H(x)-e k (t)n(x)) X k (x 1 t) 
X k (x,T) = Oi> k (x,T) 



(33) 



e k (t) = (1 - ^(t) ^Im( X k {t)\nW k (t)). (34) 
a 

Existence and uniqueness of solutions ip k and % fe of the above equations result from an easy 
adaptation of Lemma 3, as for the proof of e k ,e k € L 2 (0,T) for all k e N. 

Remark 2 : Note that this choice of optimizing sequence is not canonical. There exists 
other ways to guarantee that the condition (1) is fulfilled (see, e.g. [30]). However, this 
formulation includes many monotonic algorithms, e.g. the one by Krotov (presented in [29]) 
or by W. Zhu and H. Rabitz [34] which are often used in the numerical simulations. 

6.2 Properties of the sequence 

We present here two results about the sequence (e fc )fegN- The proofs can be found in [19, 20]. 
These results state that (e k )keN defined by (31) — (34) is bounded in L°°(0, T) and that the 
corresponding sequence (j(e fe )) feGN increases monotonically. 

Lemma 13 Given an initial field e° € L°°(0,T), let us define M by: 
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The sequences (e k )k^ and (£ fe )feeN satisfy: 

VfceN, \\e k \\ L oo {0iT) < M, \\e k \\ L ^^ T) <M. 

Lemma 14 The sequence (e k )keN defined by (31) (34) ensures the monotonic convergence 

of the cost functional J in the sense that: 



+ r - 1 lk fe+1 - 2*11 -n + 7 " 1 II** " r, (n t\ . (35) 



J{e k+1 ) - J(s k ) = (ip k+1 (T)-i; k {T)\0\ij k+1 (T)-iP k (T)) 

,2 /2_. 

Il 2 (o,t) + i, (J x y H c " c Hl 2 (o,t) 

i/iere exists Z e o swc/i iftai lim J(e k ) = l £ o. 

In order to study the convergence of (e k )kem, we will need to estimate the gradient of J at 
each point e k . Such an estimate is obtained in the next lemma. 

Lemma 15 There exists A > 0, depending on ji, O , a, 8, rj and T , such that: 

||VJ( £ fc )|| Ll(0 , T) < A (|| e fc - ? fe - 1 || L2( o,T) + II?*- 1 - e k -'\\ W) ) . (36) 
Proof: Thanks to (20), we have: 

X7J(e k )(t) = -2(ae k (t) + Im{ X k -\t)\^ k {t)) + Im( X e " (t) - X^ 1 (t)|/# fc (t)» 

= -2 (a (l - (e fe (t) - g*" 1 ^)) + Im( X e \t) - X k ~\t)\^ k (i)>) (37) 
where x £ is the solution of (16) with ?/> = tp k and e = e fe . 

Next, X £ — X fc_1 i s the solution of equation (17) corresponding to e — e k , e' = s k ~ 1 — e k and 
X = X k l - The associated estimate (18) then gives: 

\\x ek (t)-x k - 1 (t)\\v < 2|| M || L ^||e fe -^- 1 || L1(0 , T) ||O|U + ||O(^(r)-^- 1 (r))|| L2 

< A\\fi\\ L ^\\OU\\e k - t- 1 ^^ + W^- 1 - e^lUncT))- (38) 
Combining (37) and (38), we obtain (36) with A = 2y/T ^4T||0||*||/i|||oc + a (l - ^ . 



6.3 Limit points of the sequence 

We now present some result about the limit points of (e k )keN- These results give first hints 
about the relationship between these limit points and the set Cj of the critical points of the 
cost functional J. Thus, we obtain a first case of convergence. 

Lemma 16 Let (s n ) n eft be a weakly convergent sub-sequence of (e fe )fe 6 N in L 2 (0,T). Then 
(e fe ")„ e N converges in L°°(0,T) towards a critical point of the cost functional J. 

Proof: Let (e fc ")„ e N be a weakly convergent sub-sequence of (e fc )fceN in L 2 (0,T) and let 
us consider /eN. Equation (35) ensures that (e kn+e ) ne jq also converges weakly (and has the 
same limit as (s kn ) ne ^). Thanks to Lemma 5, the sequences (x fc " +£ )„ eN an d (V ;fe " +£ )„ eN 
converge strongly in the space C([0, T];L 2 ). Thus, we obtain by bilinearity the strong con- 
vergence of both sequences ((x fe "lMlV'' c "))„ eN and ((x k "\^ krl+1 )) neN in L°°(0,T). 
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According to (32) and (34), (e fe ")„ e N also reads: 



(1-S)(1- V ) e k "+u kr . 



where \v\ < land where u k (t) = -- — ^-Im (x kn (t)\M kn (t))--Im (x kn (t)\M kn+1 (t)) 



strongly converges in L°°(0,T). Note again that given fel, (ufc„+£)neN also converges in 
L°°(0, T) (towards the same limit). For all k e N, the absolute value oiuk{t) can be estimated 

^ U,(t\\<m 4 ^H^II°II* 

\u k {t)\ <m~ 



a 



Let us prove that (s kn ) n eN is Cauchy in L°°(0,T). Consider e > 0. There exists m > be 
such that 



2m£M J "<|. 
ji 

Since the sequence (wfc„-£)neN is Cauchy for all I with < I < m, we have: 
3n 2 > 0/ Vs > n 2) V<? > 0, ||ufe s+5 -^ - Mfc s -^ ||l~(o,t) < 

Let n be an integer fulfilling the conditions: 
Vp > 0, \u kn +" -u K \< 



4||£°||l~(o,t) : 



4ni ' 

k n > m, n> n 2 - 



(39) 



(40) 



(41) 



fen-l 



Let p be a positive integer. Since we have, for all n e N*, e*™ = //"e + ^ v^u^-j-i we 
obtain 



j=k„ 



711—1 



(42) 



Jl 3=0 

Thank to (39) and the two first conditions of (41): 

U^fen+P _ v k ~) £ 0\ 



\l°°(o,t) - 4' 



fen+p-l 



< 



L°°(0,T) 



j = k r . 



3 — 



According to the condition (39), the third term of (42) can be estimated by: 



fcn-i 



31 



< 2m^ < - 



L°°(0,T) 



Lastly, < 1, the third condition of (41) and the Cauchy property (40) allows us to estimate 
the last term of (42): 



ni — 1 



J=0 



L°°(0,T) 



ni — 1 
J=0 



|Ufe n+p -j-l — —J — 1 1 



L°°(0,T) 



e 

< -. 
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We have thus proved that for all e > 0, if n is large enough then, for every p > 0, 

||£ fen+p -e fen ||z,°°(o,T) < e, 

which proves that (e fc ")„ e N is Cauchy in L°°(0,T). 

We denote by e the limit of (e fc ")„ G N- Thanks to (35), (e*")„ e N also converges towards 
e. Passing through the limit in (31) — (34), we then deduce that e belongs to Cj, according 
to definition (21). 

■ 

Let us denote by C e a c Cj the set of the limit points of (e kn ) n &i- As stated in Remark 
1, for a > 6T||^||| oo ||0||*, Cj, and consequently C e o, are reduced to one point. By means of 
Lemma 13, the convergence of the sequence (e fc )feeN is then guaranteed in this case. 

Remark 3 : In addition, the uniqueness of the critical point implies that the limit in this 
case is necessarily an extremum of J. 

In order to obtain the convergence for all a > 0, we need to study more precisely the asymp- 
totic behavior of the sequence (e k )kefi in the neighborhood of C e o. A standard argument of 
compactness applied to C' £ o enables us to obtain the following result. 

Lemma 17 Let denote by d^ the distance corresponding to the £°°(0,T) norm. One has: 

doc(e fc ,C £ o)^0. (43) 

Remark 4 : By means of the monotonicity property, we find that J ~J. £ o on the set C e o 
(with l e o — linifc^+oo J(e k ) ). It is then possible to apply Lemma 12 with Cj = C\o since the 
assumption that Cj is connected is only necessary to ensure that J is constant on this set. 
It can however be proved that C e o is connected (see [25]). 



7 Convergence of the sequence 

It is now possible to prove the convergence of the sequence (e k )k£fi by a Cauchy argument. 

Theorem 2 Suppose that e° G £°°(0,T). The sequence (e fc )fcgN defined by (31) — (34) is 
convergent in L 2 (0,T). 

Proof: We still denote by l E o the value of J on C e o and by J the shifted cost functional 
J — l s o. Suppose first that Vfc G N, J(e k ) ^ 0. By (43), there exists k such that (30) holds 
(with Cj = C £ o) for all e k with k > k . Consider an integer k > fc . We have: 

((J(e k )f - (J(e k +^ > I (J(e k ^)-J{e k )) (44) 



> "A (f 2 lW+1 F ~fc||2 



+ (£_l)||g*_ e fc||2 2(0iT) ) (45 ) 



> ^^(l| £ fc+1 -e fe |U 2( o,T)+l|e fc -e fe || i2(0 , T) ) (46) 



> 



K6a (S,r,) | U fc+l _ k 



e — e l|i 2 (0,T), 
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where a:g v \ = — r- — — . The inequality (44) comes from the concavity of s i— » s , 

' ; max(d, 77) 2 

whereas (45) is a consequence of (30) and (35). Inequality (46) follows from (36). 

Since | (^J(e k )j j is a Cauchy sequence (as a monotonic sequence bounded by (2||O||») ), 

we obtain that (e fe )fceN is also a Cauchy sequence. 

If there exists k\ such that J(e kl ) = 0, the monotonicity of the algorithm implies that 

J(e kl ) = J(e kl+1 ) = J(e kl+2 ) = ... 
and by (35) the sequence (e fc )fc C N is constant for k>k\. 



Remark 5 : Thanks to the definition of the sequence (e k )kq_n and to the regularity of the 
solutions ^ andx of the appropriate Schrddinger equations (see lemmas 1 to 3), we can easily 
prove by induction that if e° G W hoc (0,T), then for all k E N, e k e W 1,oo (0, T). 



8 Rate of convergence 

The rate of convergence can be now evaluated by a second use of the Lojasiewicz inequality. 
The result is summarized in the next theorem. 

Theorem 3 Let us denote by e°° , the limit of (e fc )fc 6 N defined by (31) — (32) and 9, k the 
real numbers appearing in (30), where C e a — {s°°}. 

~ 1 , L_ 

If 6 < —, then there exists c > such that \\e k — e°°||L 2 (o,T) < c k 1 - 2e ■ 

If 9 = —, then there exist c' and t such that: 
2 

l|£ - £ IIl 2 (o,t) < c e . (47) 
Proof: As in the proof of Theorem 2, let be ko, an integer such that 

W>fc ||VJ( £ £ )llL 1( o,T)>«|J(^)r^ (48) 
Let us fix k > ko and introduce A k defined by: 

oc 

e=k 

With no loss of generality, we may assume that A fc > for all k > ko. Summing (46) between 
k and +00, we obtain: 

A 

This estimate, combined with (48), with £ = k yields: 



l|VJ(e fe )|| Ll(0 , T) >«(^^A fc )" 



From Lemma 15, we obtain: 



A(A fc - 1 -A fc )>«(^^A fc ) 1 ", 
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which may be written as follows: 

A fe-i _ A k 

19 ~~ x ^ ~ 1 

with (3 = — = — and v = j ( K a ^'" } J 6 . Suppose now that 9 = -,i.e., (3=1. The equation 

(49) then becomes: 

(l + v) k °A ko (—!—) k > A k , 

and (47) is proved with c' = (1 + v) k ° A k ° and r = ln(l + v). 
Suppose now that 9 < \. Let be r e]0, 1[, and suppose first that: 

( A fc)/3 > r (A fc ~ 1 )' 3 . 
Since 1 — (3 < 0, the function s s 1 ~' 3 is concave and we have: 

Afe-l _ Afc A^^ 1 — A fc 

{Ak y-P _ ^k-if-fl > (/3 _ 1) (Afc _ i)/3 > (/3 - l)r ^ >[fi- l)rv. 
In the other case: 

( A fe)i-/3 _ (A fc - 1 ) 1 -' 3 > (A fc ) 1 "' 3 - (r^A*) 1 "' 3 = (1 -r^XA*) 1 -' 3 > (1 - r i ^)(A' C0 ) 1_/3 . 
Thus, in any case, there exists v' > independent of fc, such that: 

(A fe ) 1-/3 - (A fe_1 ) 1_/5 > w'. (50) 
Consider now k' > k, the inequality (50) implies that for a small enough c, one have: 

— i 

A k ' < (y'{k' -k) + (A k ) 2 -^ l ~ 2S < ck'-^Z, 
and the result follows. 



Remark 6 : Thanks to Lemma 10, we have thus obtained that if a > 6T||^|||°o||0||* the 
convergence of the sequence is at least linear. 
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Appendix: Proof of Lemma 11 

Consider e e L 2 (0,T) and J defined by (14). For reason of simplicity, we suppose that 
J(e) = 0, V J(e) = 0. 

Thanks to Lemma 9, the operator Hj(e) is a Fredholm operator. The Fredholm alternative 
states then that either Hj(s) is bijective or KerHj(e) = span(ipi, ...,ip rn ), with m > 0. Let 
us denote by n, the orthogonal projection on KerHj(s) (with n = if KerHj{e) = 0). The 
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operator L = U + Hj(e) is then bijective on L 2 (0, T). 

We are now in the position to apply the local inverse mapping theorem to C = II + V J 
(analytic version, see [33], Corollary 4.37, p. 172), that asserts there exist V and V two 
neighborhoods of in L 2 (Q,T) and K : V' — > V an analytic mapping such that: 

Ve' g V, K(C{e')) = e', Ve" e V, C(K{e")) = e" . 

Since C and K are C°°, there exist C and C such that: 

Vei,e 2 eV, ||£(£ 2 ) - A £ i)IU 2 (o,t) < C||£ 2 - EiIIl^o.t) 
Vei.e^e^', ||^(4) -A-(ei)|| i2(0 ,D < C"II4 - 4 IU-(o^)- 

Consider now e'eFnf For ( e M m such that Y!?=i Cjfj ^ v > lct us dcfmc r : C ^ 
J{K(J2™ =1 (j<Pj)), and £ e M m such that lie' = EjLi^^- Let us first estimate Vr(£). 
Using lie' e U', we obtain: 

|vr(OI < c"||VJ(^(n £ '))IU 2 (o,T) = c"||vj( £ ') + vj(^(n £ 0)-vj( £ ')ll^(o > T) 

< C"(\\VJ(e')\\ LH0 , T) + C\\K(Ue') - e'\\ L 2^ T) ) 

= C"{\\VJ(e')\\ L 2^ T) +C\\K(Ue')-K(Ue' + VJ(e'))\\ L 2^ T) ) 

< c||VJ(e')ll^(o,T), 



(51) 



where c = C"(l + CC). On the other hand, one has: 



|j( £ ')-r(OI - |J( £ ')- ^(^(n e '))l 



Z* 1 ijf £ ' + s (]f(n £ ')-£'))^ 
J Q as \ ) 

y (Vj(e' + s(X(ne') - e')) , K(ILe') - e') ds 

< \\K(ne>)-e')\\ L 2 (0 . T) f \\VJ(e')\\ L 2 (0 . T) +Cs\\K(ne>)-e>\\ L 2 (0 . T) ds 

Jo 

= \\K(ILe') - e')\\ W) (\\VJ{e')\\ W) + ^\\K{Iie') e'\\ W) ) 

< C||VJ( e ')lli a( o,r). 



cc' 



(52) 



). By diminishing V, the Lojasiewicz inequality (29) applied to the 



where d = C'(l 

analytic functional T states that there exist 9 G]0, 1/2], a > such that: 

|vr(OI > TOI 1-9 = \J{e') - r (£) - AOI 1 "" 



> 2l J ( £ ')| 1_ 



-Ij^O-r^)! 1 - 9 . 



Combining (51) and (52), we obtain: 



C ||vj( £ ')IIlw) > 2 |J(£ ' )|] 



c'\\VJ(s'^ 2{1 - e) 



1L 2 (0,T)' 



and the result follows. 
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